Spectral functions of strongly interacting isospin-| bosons in one dimension 
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We study a system of one-dimensional (iso)spin-i bosons in the regime of strong repulsive inter- 
actions. We argue that the low-energy spectrum of the system consists of acoustic density waves 
and the spin excitations described by an effective ferromagnetic spin chain with a small exchange 
constant J. We use this description to compute the dynamic spin structure factor and the spectral 
functions of the system. 
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Physics of one-dimensional Fermi systems has long at- 
tracted the interest of both theorists and experimental- 
ists. Interactions between particles have a strong effect 
on the properties of these systems. Interacting fermions 
form the so-called Luttinger- liquid state [l[, whose ex- 
citations are bosons with acoustic spectrum, e(q) oc \q\. 
Recently it has become possible to confine ultracold gases 
of bosons to elongated traps [2, S|, effectively creating 
systems of one-dimensional bosons. The properties of in- 
teracting one-dimensional spinless bosons are in many re- 
spects similar to those of spinless fermions. In particular, 
they too form a Luttinger-liquid state at low energies. 

In a recent experiment [i| bosons with two internal 
degrees of freedom, which can be viewed as compo- 
nents of (iso)spin-i, were confined to one dimension. 
For spin-i particles the difference between the Bose and 
Fermi statistics is of fundamental importance. Indeed, 
spin-independent interactions between one-dimensional 
bosons favor ferromagnetic spin ordering_l5| , whereas for 
fermions the ground state spin is zero [6(. As a result, 
the low-energy spin excitations of the boson system are 
magnons with quadratic spectrum e{q) oc q 2 , and the 
system is no longer a Luttinger liquid. 

In the absence of the effective theory of interacting 
spin-i bosons in one dimension, considerable progress 
has been made recently by focusing on the regime of 
very strong repulsive interactions 0, Q . In this paper 
we show that this regime allows for a remarkably simple 
theoretical description, in which there are two types of 
low energy excitations: acoustic density waves and the 
spin excitations described by a one-dimensional Heisen- 
berg model with a very small ferromagnetic exchange 
constant J. The theory is applied to the calculation of 
the dynamic spin structure factor and the spectral func- 
tions of the system. Unlike Refs.0 andd, our conclusions 
are not limited to spin excitations of small momentum 
q — > 0. In addition, although the frequency lu is assumed 
to be small compared to the typical kinetic energy of the 
bosons Ep ~ (Tin) 2 /m, it can be of order of the small 
exchange constant J. (Here n is the one-dimensional 
density of bosons and m is their mass.) 

The model we consider is that of one-dimensional 



(iso)spin-i bosons interacting with repulsive spin- 
independent potential V(x — y). For simplicity, we con- 
centrate on the most realistic regime of short-range inter- 
actions, V(x — y) = gS(x — y); the generalization to the 
case of finite-range repulsion is relatively straightforward. 
The strong repulsion regime is achieved at 7 3> 1, where 
7 = mg/H 2 n is the dimensionless interaction strength. 

As the first step, we show that at low energies the ex- 
citation spectrum of the system consists of independent 
phonon and magnon excitations. This effect is essen- 
tially equivalent to the well-known spin-charge separa- 
tion in interacting one-dimensional electron systems [Io| . 
Our arguments follow the discussion ll|, [l2j of that phe- 
nomenon in the limit of strong repulsion. 

In the Tonks-Girardeau limit 7 — > +00 the repul- 
sion effectively forbids any two particles to occupy the 
same point in space, regardless of their spin. Thus 
the density excitations of the system are those of a 
gas of spinless hard-core bosons, or, equivalently, those 
of non-interacting gas of spinless fermions 0, where 
the same constraint is enforced by the Pauli principle. 
It is convenient to treat the low-energy excitations of 
one-dimensional spinless Bose and Fermi systems in the 
framework of the hydrodynamic approach [l[ and write 
the Hamiltonian in the form 



Hph — 



[K{d x 9f + K- 1 (d x <t ) f}dx. (1) 



Here cj> and 9 are bosonic fields satisfying the standard 
commutation relation \4>{x), d y 9(y)] = iit8(x — y). The 
Luttinger liquid parameter K and the phonon velocity u p 
are determined by the interactions. In the case of hard- 
core bosons K = 1, while the effective "Fermi velocity" 
Up = irhn/m. 

In the limit 7 — > +00 any collision of two bosons results 
in perfect backscattering. As a result the bosons become 
distinguishable particles. Indeed, if boson 1 is to the 
left of boson 2, i.e., x\ < xi at some moment in time, 
then this property cannot be changed as a result of any 
collisions between particles. Thus one can number all 
particles by an integer I in accordance with their positions 
along the cc-axis. In this limit the spins of the bosons 



2 



do not interact, and each state of N bosons is 2 Ar -fold 
degenerate. A coupling of the spins appears only when 7 
is finite. At 7 ^> 1 a collision of two bosons, I and I + 1, 
may result in their forward scattering, in which case the 
particles exchange their spins. Since for spin-^ particles 
the spin permutation operator P/,z+i can be expressed as 
Pi,i+i = 2Si ■ Si + i + 1/2, this gives rise to coupling of the 
spins of the nearest-neighbor particles: 



H„ 



1 



JSi-S, 



i+i- 



(2) 



Thus at 7 ^> 1 the low-energy excitations of the sys- 
tem are given by the acoustic phonons, described by the 
Hamiltonian |T]) and the spin excitations of the Heisen- 
berg Hamiltonian ((2|). 

A similar separation of the density and spin excita- 
tions is well known in the case of strongly interacting 
one-dimensional fermions, where it was first derived [131 ] 
from the exact solution of the infinite-!/ Hubbard model. 
The sign of the exchange constant J is determined by the 
requirement to either symmetrize or antisymmetrize the 
wave function with respect to the permutation xi <-> a;j+i; 
the coupling is antiferromagnetic for fermions, J < 0, and 
ferromagnetic for bosons, J > 0. On the other hand, the 
magnitude of the exchange constant J is determined by 
the amplitude of the forward scattering of two neighbor- 
ing particles, regardless of their statistics. Thus we find 
the same value of J as in the case of fermions with strong 
short-range repulsion, 



J 



2tt 2 h 2 n 2 
3 rwy 



(3) 



see Eq. (22) of Ref. [TT|- The effective theory (0), © of 
the spin subsystem is consistent with the recent thermo- 
dynamic Bethe ansatz results 14 j. 

The Hamiltonian describing all the low-energy excita- 
tions of the system is the sum H p h + H a . An important 
assumption in its derivation was that all the relevant en- 
ergy scales in the problem, such as the temperature T, 
are small compared to the bandwidth (the "Debye fre- 
quency") of the phonons Ep . In the following we limit 
our discussion to the most interesting case of T = 0. 

The ground state of the ferromagnetic spin chain ([2]) is 
fully spin-polarized. The excitations near this state, the 
magnons, have the well-known spectrum 



e(Q) = J(l-cosQ), 



(4) 



where Q is the wave vector defined with respect to the 
lattice of the spin chain ([2]) and varying in the range 
— 7r < Q < tt. Since the spins are attached to par- 
ticles filling the real space with density n, the phys- 
ical momentum of the magnon is p = HnQ [TR fig . 
In the limit of small p, the spectrum ([?]) is quadratic, 
e (p) = p 2 /2m*. Using Eq. ([3]), one finds the effective 



mass m* = (3/27r 2 )7m, in agreement with the result of 



Ref. |17|. 

Let us now illustrate our approach based on the sepa- 
ration of the density and spin excitations in the form (|TJ) 
and ([2]) by calculating the dynamic spin structure factor 



S±(q,u) 



dx dt 

"ST 



,— iqx-\~iojt j o-f 



(S+ (x,t)S- (0,0)). (5) 



Here S(x) is the spin density operator, S ± — S x ± iS v , 
and the expectation value (. . .) is evaluated in the fully 
polarized ground state of the system, with the polariza- 
tion assumed to be directed in the positive z-dircction. 

We start by expressing the spin density operator S(x) 
in terms of the particle density operator n(x) and the 
spin operator Si, 



S(x) = n(x)S, 



l(x)- 



(6) 



Here l{x) is the operator of the number of particles to 
the left of point x, i.e., d x l{x) = n(x). Its presence in 
Eq. ([6]) accounts for the fact that the operator S(x) acts 
on site I of the spin chain ((2|) attached to the boson at 
point x, cf. [l2|. 

The problem of zero-temperature properties of 
strongly interacting bosons is considerably simpler than 
that of fermions [l2j, because of the simplicity of the 
ground state of the ferromagnetic Heisenberg model ([2]) 
and its single particle excitation spectrum ((4|). In par- 
ticular, the correlator (S' ; + 5 ; 7) C r for the spin chain ([2|) is 
easily found as 

dQ iQ(l-l')-iQ(Q)t 



(5+(t)5,7(0)) ff 



2vr 



(7) 



where ft(Q) = e{Q)/Ti is given by Eq. g]). Substituting 
Eq. ([6]) into ([5]) and using ([7]), we find 



dx dt dQ 



g — iqx+i[ul — fl(Q)]t 



x (e^^-^nix^MO^)) . (8) 

The expectation value (. . .) p h is performed in the ground 
state of the phonon Hamiltonian ([1]). To evaluate it, we 
use the standard hydrodynamic expression for particle 
density n(x) = n+ ^d x <fi(x), and the resulting expression 
for the particle number 



l(x) — nx H — 4>{x). 



(9) 



In the low-energy limit one can neglect the d x 4> correction 
to n(x), and replace it with the average value n. However, 
it is important to include the field <j) in Eq. ([9]) when 
evaluating the exponential in the second line of Eq. ([8]) . 
The latter calculation is performed using the standard 
techniques [l(, resulting in 



3 iQ[l(x,t)-l(0,0)}\ — 



a inQx 



'ph [(1 + iDt) 2 + (Dx/u p ) 2 ] (Q/27T)2 



(10) 
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where D ~ Ep/h is the phonon bandwidth. 

In the denominator of Eq. IjlOp one can neglect x com- 
pared to u p t. Indeed, to this approximation one finds 
that Eq. (HUD falls off at Q ~ 1/y/MDt), resulting in the 
estimate x ~ 1/nQ ~ ^\ii(Dt)/n < u p t, cf. [IIH- The 
remaining calculation is straightforward, and one finds 



- n{q/n)) 



lu — Q(q/ri) 



D 



q 2 /(2K 2 n 2 )-l 



T{q 2 /2ir 2 n 2 ) 1 

(11) 

Here is the unit step function. Its presence in 

Eq. (jlip expresses the obvious fact that the minimum 
energy of a spin excitations with momentum q is e(q/n), 
Eq. Q. 

The structure factor is essentially a Fourier trans- 
form of the correlation function Gj_(x,t) discussed re- 
cently by Zvonarev et al. 0. Their treatment is lim- 
ited to the regime q <C n; in which case our results are 
consistent with Eqs. (13) and (14) of Ref. |8|. On the 
other hand, our calculations show interesting behavior 
at larger q, especially the additional features at uj <C J/h 
and q ~ ±27m, ±47m, . . .. 

We now apply our technique to the calculation of the 
single-particle spectral functions of the system 



dx dt 
2tt 

dx dt 
2tt 



e -M^ s ( M )^t( 0) 0)), (12) 
e-^ +iwt (^(0,0)^(x,t)), (13) 



where 'i/'s (x) is the annihilation operator of bosons with 
spin s = f,|. 

As discussed above, at strong repulsion (7 — > +00) 
the density excitations of the system are identical to 
those of a gas of spinless hard-core bosons, whose density 
^{x)^{x) equals the true particle density n(x). (Here 
if? is the annihilation operator of the hard-core bosons.) 
Then, assuming that the ground state is polarized in the 
positive z-direction, one concludes that operator ip^ sim- 
ply destroys a hard-core boson, i.e., ip-\{x) — if?{x). In 
the low-frequency regime oj <^ D the spectral functions 
(q,uj) can then be obtained in the framework of the 
hydrodynamic approach based upon the Hamiltonian ((T|) 
with K = 1. In this method the annihilation operator 
if? is expressed in terms of the bosonic fields entering the 
Hamiltonian (JTJ) as 



i6{x) y^^ & i2-Kjl{x) 
3=1 



-i2Trjl(x) 



(14) 

Here one should use the hydrodynamic form §§§ of the 
particle number operator l(x). 

Compared to the first term in the right-hand side of 
Eq. (|14| . the remaining ones are formally irrelevant, i.e., 
their contribution to the observable quantities is expected 
to show additional power-law suppression at low ener- 
gies. The reason for writing the full expression (| 14[) is 



that this form accounts for the discreteness of particles 
by enforcing the condition of l{x) being integer [l|, l22l |. 
As a result, at lu <C D ~ nu p the spectral function 
A^(q,oj) = A^(q 7 uj) + A^(q,u>) shows not only the ex- 
pected feature near q — 0, but also weaker features at 
q = ±27m, ±47m, . . ., 



E 



pooA\j\ 0(uj 2 - u 2 (q - 2irjn) 2 ) 



r((j 



^) 2 )r((j + ^) 2 ) 



\oj — u(q — 2irjn)\ 
2ixnu 

\uj + u(q — 2irjn)\ 
2-nnu 



U+i 



(15) 



The hydrodynamic approach does not enable one to ob- 
tain the numerical coefficients and Aj. To find 
them, one can compare the equal-time Green's function 
com put ed within this approach with the exact results 
[3 lH EH- This results in Poc = 0.92418, A = 1, 
Ai = l/16, A 2 = 9/2 16 ,.... 

In this paper we are primarily interested in the spectral 
function A^ , because unlike A^ , it is sensitive to the non- 
trivial spin properties of the system. (The other spin- J, 
spectral function, Aj, obviously vanishes.) To evalu- 
ate A~^(q, oj), one needs to express the operator ipi(x) 
in terms of the density and spin variables entering the 
Hamiltonians {1} and ([2]). Following the ideas of Refs.fTil 
and we identify 



i'lix) = ^(x)Z l{x)d . 



(16) 



The presence of the hard-core boson operator 'J accounts 
for the change in the total number of particles in the sys- 
tem, when a particle with spin-| is destroyed. In addi- 
tion, the number of sites in the spin chain ^ reduces 
by one. This effect is accounted for by the operator Zi d , 
which by definition removes a site at position I in the 
spin chain, provided that the spin at that site is J,. (Oth- 
erwise, the result is zero.) 

In the fully spin-polarized state of an infinite spin chain 
@, the correlator (ZiiZj,,)^ coincides with the spin-spin 
correlator 0. Then the substitution of Eq. (16]) into (fl~2"l) 
gives 



A+(q,Cj) 



dx dt dQ 



-iqx-\-i[oj — Q,(Q)]t 



x (^(x,t)e tQ[l( - x ^- l( -°' a)] ^(0,0)) .(17) 

\ / ph 



To evaluate the expectation value in the ground state 
of the Hamiltonian we use the hydrodynamic the- 
ory expression (fT4]) for the hard-core boson operator. 
Upon substitution of Eq. (fl4| into (fP7|) . the effect of the 
j ^ terms amounts to the extension of the range of Q- 
integration from (— tt, 7r) to (—00, +00). Then the corre- 
lator in the second line of Eq. (fl~7|) is computed with the 



4 



help of the relation 



-i6(x,t)+iQl(x,t) i6(p,0)-iQl(0,0)\ _ 

I ph 



[iD(t - x/u p ) + 1] A Q [iD(t + x/u p ) + 1} X Q 



(18) 



with Aq = (Q/tt± l) 2 /4, obtained using the standard 
techniques [l|. 

Similarly to our derivation of the dynamic spin struc- 
ture factor (fTTj) , at low frequencies one can neglect the 
x-dependence in the denominator of Eq. (|18[) and find 

<&(u -Cl(q/n)) 



T(q 2 /2TT 2 n 2 + 1/2) 



D 



to — Q(q/n) 



D 



(j 2 /(27T 2 n 2 )-l/2 



(19) 



This expression is the main result of our paper. It 
is worth noting, that similarly to the case of Ay(q,u), 
Eq. (|15[) . the hydrodynamic approach does not enable 
one to accurately determine the prefactor in Eq. (fill)) , 
whose calculation at this time remains an open problem. 

The spectral function A~^(q, to) is defined as the Fourier 
transform (Tj"2"]) of the spin-| boson Green's function. 
The latter was discussed recently by Akhanjee and 
Tserkovnyak Their theory focused on the Jt — > 

oo limit, and accounted only for the long-wavelength 
magnons, q -C n. Calculating the inverse Fourier trans- 
form of Eq. (|19[) under these assumptions, we get 

, + , n 1 ( ihn 2 x 2 " 



y/2nDJ/h 



1 

it + 



exp 



2Jt 



. (20) 

Comparison with the considerably more complicated 
Green's function given by Eq. (7) of Ref. shows the 
same oscillating exponential factor (up to a missing it 2 
in their exponent). Further, assuming \x\ <C u p t in the 
result of Ref. 0, we find that their prefactor is consistent 
with our Eq. (|20| . 

It is interesting to compare the spectral function (|19p 
with that of strongly interacting fermions [12] • The latter 
calculation, performed in the limit J <C hu), shows the 
same Gaussian peak as a function of q at small u> as the 
expression (|19p at fl cx J — > 0. In both cases the peak 
gives the leading contribution to the density of states, 
obtained as g-integral of Af(q,uj), resulting in v(co) 



\/\J(jJ hx{p/u)), cf. [13, [23j] . In addition to the peak at 
q = 0, the spectral function of the fermion system shows 
weaker features at the Fermi surface, q — ±kp, as well 
as the shadow-band features at ±3fcp, ±5fcp, etc., with 
the Fermi momentum kp — irn/2. At J/lo — > the 
boson spectral function ([T9")) does not show any additional 
features. However, at lo < J/h we find a sequence of 
additional features at q = ±2nn, ±4-7m, etc. 

To summarize, we have developed a new approach 
to study the low-energy properties of a gas of one- 
dimensional (iso)spin-i bosons with strong short-range 



repulsion. Our method is based on the separation of den- 
sity and spin variables in the form (|TJ) and ([2]) and the 
expression (fT6)) for the boson annihilation operator. We 
applied this technique to the calculation of the dynamic 
spin structure factor (fTTj) and the spectral function (|19[) . 
At small lo they both show Gaussian peaks as a function 
of q centered at q — 0, as well as sequences of additional 
features at q = ±2nn, ±47m, . . .. 
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